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FOREWORD 


This, Volume I of the report "Spherical Roller Bearing 
Analysis" documents the analysis, program design and algorithm 
detail employed in the generation of the computerized analytic 
design tool SPHERBEAN. Efforts involved in the generation of 
the code were accomplished under contract NAS3-20824 issued by 
the NASA-Lewis Research Center of Cleveland, Ohio and admin- 
istered by the Structures and Mechanical Technologies Division. 
Funding was provided by the Product Assurance office of the Army 
Aviation Research and Development Command, St. Louis, Missouri. 
The technical monitor was Mr. H. H. Coe. The work was performed 
at SKF Industries, Inc., King of Prussia, Pennsylvania, during 
the period June 1978 to December 1980. 

Technical project leadership was executed by Mr. R. J. 
Kleckner, with contributions from Dr. J. Pirvics and Messrs. 

G. Dyba and T. Deromedi. 
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NOMENCLATURE 

C Cage rail/ring land radial clearance (M) 

6 

Cjj Drag coefficient for fluid flow over a cylinder 

E' ( 1-Vl + 1 ~ v i ) , (N/M 2 ) 

Er 7T E 2 

E^,E 2 Elastic moduli for bodies 1 and 2 (N/M^) 

^A Lubricant parameter given in [19] 

G Lubricant coefficient of thermal expansion (l/C?) 

h 0 Film thickness (M) 

AAA 

i,j,k Unit vectors along the x,y,z axis, respectively 

K l Lubricant parameter given in [18] 

M Roller mass (kg) 

n Number of rollers per row 

n Number of cage guide rails 

o 

n r Number of rows of rollers 

n g Number of slices 

n^ Total number of rollers 

q Heat generation rate (watts) 

r Equivalent cage rail radius (M) , see equation (25) 

6 

V,U Sliding and entrainment velocity vectors (M/sec) 

w Equivalent width of cage guide rail (M) , see equation 

8 (26) 

a Pressure-viscosity coefficient at elevated temperature 

(M 2 /N) 
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NOMENCLATURE (continued) 


e 

Eccentricity 


n 

0 

Lubricant viscosity at ambient pressure, and 
temperature (N-sec/MO 

elevated 

V 

Kinematic viscosity (CS) 


V , V 
1 2 

Poisson’s ratios for bodies 1 and 2 


p 

Lubricant density at elevated temperature (gm/cm) 

P' 

P 

1 6 

Curvature sum, ^ ^ (1/M) 

r l 2 

Lubricant density at 16°C (60°F) , (gm/cm 5 ) 


a 

Contact stress (N/M^) 


T 

2 

Lubricant shear stress (N/M ) 


T c 

Critical shear stress, 6.89 x 10^ N/M^ (1000 

psi) 

*s 

High contact stress factor 


SUBSCRIPTS/ SUPERSCRIPTS 


i 

Refers to the i-th roller, i = 1,2 

n t 

j 

Refers to the j-th row of rollers 


k 

Refers to the k-th slice, k = 1,2,...., 20 


m 

Refers to the m-th raceway, m = 1, outer, m 
inner. 

= 2, 

NOTATION 



a 

A vector with components a , a , a 

x y z 


1 a 1 

Magnitude of vector a 


(x) 

Vector cross product, a x b = c 


C) 

Vector dot product, a • b = d 
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1.0 INTRODUCTION 

The general class of load support systems (LSS) consists 
of mechanical assemblies which transfer force and moment vec- 
tors. The transfer is constrained to maintain component 
displacements within prescribed, and frequently severe limits. 

Shaft mounted rolling element bearings form a subset of 
the general LSS class. They are characterized by high load 
carrying capacity, low power loss, relative insensitivity 
to load fluctuation, rotation reversal and tolerance for start 
stop operation and reversal of rotation direction. 

Various rolling element bearing configurations have 
evolved to service specific application requirements within 
this subset. The angular contact ball bearing, for example, 
supports combined loading and tolerates some misalignment. 

Because rolling elements with "line" as compared to "point" 
contact conjunctions offer superior capacity for a given design 
volume, bearings having tapered and cylindrical roller geom- 
etries have evolved to support large loads. 

The relative displacement of components is a physical reality 
wherever mechanical assemblies transfer force. Structural 
elements distort within the generally asymmetric assembly and 
cause the misalignment of bearing raceways. This departure from 
the idealized "rigid" assembly can be controlled but not 
eliminated by manufacture, assembly and operation. In the 
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presence of such deflections, and the associated misalignment, 
line contact geometries can be optimized only for a specific single 
load condition [1]*. Severe penalties, in bearing life and oper- 
ating performance, are sustained during less than optimum load 
support. 

The self aligning spherical rolling element bearing, Figure 1, 
answers some shortcomings of the bearing types noted above. The 
geometry is unique. At low loads, the load vector is transferred 
by point contacts. At higher loads, modified line contacts per- 
form this function. The bearing also supports combined radial 
and axial loading. This versatility has led to successful imple- 
mentation in the large load support systems required by steel , 
paper and marine industries. Smaller mechanical assemblies, such 
as planetary gear reduction sets, have also seen successful 
application. 

The rolling element bearing subset of the LSS class, in 
today's technology, is being required to operate at ever increas- 
ing DN values. Ball bearings, for example, have seen numerous 
high-speed applications up to the 2 X 10^ DN level. Cylindrical 
and tapered roller configurations are being asked to follow in 
this regime. These new demands result from the requirements posed 
by advrnced hardware missions and the increased emphasis on 
extracting maximum energy from a given process cycle. Basic 

* Numbers in brackets denote references at end of report. 

2 

SKF TECHNOLOGY SERVICES 

SKF INDUSTRIES. INC 



41 BSC -1 -R100 


ATS IDO 06 


thermodynamics, particularly in mobile power plant design, 
dictates higher temperatures, stresses and speeds. Simultan- 
eously, the assembly is required to occupy a decreased volume 
and weigh less. The combination of these parameters defines 
a lighter assembly, under increased str ;s, in a high temp- 
erature environment. Bearings, residing in assemblies of 
decreased structural rigidity, must sustain high, combined 
loads under conditions of misalignment and furthermore, do so 
at higher speeds. 

The conventional spherical rolling element bearing design 
meets all but one of the challenges posed by these emerging 
requirements. Operating speed has been restricted to maximum 
values on the erder of 5000 rpm. DN values have peaked at 
about .25 X 10^. Efforts are now under way to reach higher 
speeds. Particular emphasis is placed on reaching a 1.0 X 10 6 
DN value. 

The ne^i to extend the operating DN regime for spherical 
roller bearings ir. new as well as current applications re- 
quires a realistic assessment of current methods for their 
design. Examination reveals that design relies on "rules", 
hand calculations, and some modest computerized simulations. 
The presence of mandated safety factors in rules for success- 
ful application reveals the measure of design performance 
buffers. At the same time, it reveals an opportunity for 
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increasing performance and thus DN values to satisfy emerging 
load support needs. 

Practical use of design reserves requires a more detailed 
understanding of, and the ability to predict, bearing perfor- 
mance within a load support system. The complexity of the 
interactions between the LSS and its environment requires 
an analytic/design tool to accurately describe the thermomechan- 
ical dialogue present [2]. Such a simulation tool can be 
created in the form of high speed digital computer software. 

Several investigators have addressed the analysis of 
spherical rolling element bearings thus setting the stage for 
the required computerized simulation. Recently, Kellstrom [3] 
explored the fundamental mechanics which control symmetric 
roller behavior. He specifically addressed roller "self 
guidance" and explored the optimization of their skewing 
angles to minimize heat generation. Wieland and Poesl [4] 
have presented an interpretation of the empirical state of the 
art in spherical roller bearing design and application. Harris 
and Broschard [51 as well as Liu and Chiu [6] have examined 
these bearings in planetary gear applications in earlier 
investigations. Palmgren [7 ] and Harris [8 ] touch on computa- 
tional procedures. Manufacturer catalog data and popularized 
applications articles serve to further highlight the need for 
a thorough analytic examination of the coupled phenomena which 
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occur during spherical roller bearing operation. 

This report documents the mathematical foundation for 
the analytic/design tool SI rtERBEAN (SPHERical BEaring ANalysis). 
Specific algorithms, deviating from those of its cylindrical 
bearing sibling CYBEAN [9] and parent SHABERTH (SHAft BEaRing 
THermal analysis) [10], are described. Program architecture 
and solution methodology are likewise described which lead to 
SPHERBEAN'S capability to address quasidynamic equilibrium 
with consideration of EHD and HD forces at the raceway and 
flange, cage skew control, heat generation and centrifugal 
loads for single or double row designs having up to 30 rollers 
per row. 

In Volume II [11] of this report, examples of program 
use are displayed with a description of computer resource 
needs. In Volume III [12] spherical bearing performance 
predicted by SPHERBEAN is compared to data obtained from full 
scale hardware tests. 


5 
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2.0 COMPUTATION OF BEARING PERFORMANCE AT CONSTANT TEMPERATURE 

Assume that a spherical roller bea~ ; g is to be simulated 
and that its performance is to be predicted as a function 
of physical dimensions, material properties and operating con- 
ditions. Then, the problem addressed is that of generating a 
computerized tool, specifically a tool which takes the form 
of a mathematical analog to a generalized physical config- 
uration. 

Consider a lubricated spherical roller bearing with its 
outer ring rigidly mounted in a supporting structure . The 
inner ring accommodates a three dimensional load vector P im- 
posed upon it by a shaft spinning at constant speed u. Lubri- 
cant is applied at a known rate and occupies a fraction, s, of 
the free space in the bearing cavity. At a given instant in 
time, the bearing components and lubricant achieve the temper- 
atures T c (°C), c = 1,2 ,8. Here the subscript c refers to 

a specific component: 

c = 1, outer ring c * 5, flange, row 2 

c = 2, inner ring c * 6, bulk lubricant 

c * 3, rollers c ~ 7, shaft 

c = 4, flange, row 1 c * 8, housing 

The set of 13 quantities P, w, C, and T c serve to define 

J — 

The planetary gear application is addressed in Section 2.4. 

6 
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the operating condition of the bearing. 

Their values are considered known and sufficient to completely 
describe the bearing application. 

Bearing response to a particular set of operating condi- 
tions is defined in terms of relative displacements and 
speeds. 

An arbitrarily chosen reference state, from which component 
displacements are measured under applied loading, is shown in 
Figure 2. It is assumed that the outer ring, and therefore the 
outer ring reference frame, is fixed in space. The cage, inner 
ring and outer ring coordinate frames are coincident in this 
reference position, and the rollers are equally spaced and 
radially positioned so that at the point of closest approach 
to either raceway, a clearance equal to 1/4 the diametral 
clearance exists. Individual rollers are identified bv their 
corresponding azimuthal position, 4. 

Inner ring response, shown in Figure 3, consists of pure 
translation, and is defined by the vector a. The consideration of ring 
rotation has been omitted because of the self aligning capability of the bearing. 

Roller displacement response is described by translation 
in a plane along the vector b, and rotations about two axis 
defined by the quantities y g (skew angle) and Y t (tilt angle), 
Figure 4. In the current work, values for the displacement 
terms b , b , y , y , with the roller speed vectors w and a; 

X U * t T — o 
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(Figures 3 and 4) specify the roller response. 

Cage response is described by a rotation about its x-axis, 
i|>, aid translation in its radial plane, c tJt c z . Figure 5. If 
the rearing contains a split (two piece) cage. Figure 6, each 
half of the cage is permitted three degrees of freedom. The 
vecor 

X * (b^, b^, Y s , Y^» |w^|, c ^ ,c 2 * a z 

i = 1,2, ,n t 

is the vector of unknowns and contains the bearing independent 
variables. Their values describe the bearing’s response to a 
given set of operating conditions. 

Values for the independent variables are obtained as 
follows. First, field equations of force and moment static equi- 
librium are formulated for bearing components (rollers, cage and 
inner ring) in terms of their displacements and speeds. Con- 
stant quantities in the field equation set characterize bearing 
geometry, material properties and specific operating conditions. 

Because equations of static equilibrium are applied in 
the presence of acceleration, quasidynamic forces and moments 
(see page .5) are included to account for the effects of such 
transient loading. 

A modified Newton -Raph son iterative numerical method is 
used to obtain the unique set of solution values for the independent variables. 

8 
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These values in turn enable the computation of derived quantities, 
such as heat generation rate and fatigue life. 

2.1 Preliminary Computations 

A substantial amount of computation time can be saved if 
those terms which are not functions of the independent variables 
are evaluated prior to entering iterative solution loops. 

Terms, such as pitch diameter, are therefore computed once and 
then stored for later use. Computation details having practical 
importance in the analysis are discussed below. 

Changes in shaft/inner- ring and outer- ring/housing fit 
pressures are computed as a function of the initial inter- 
ference, operating temperatures and speed. The changes in 
diametral clearance, from off-the-shelf to operating value, 
are evaluated as functions of operating fit pressures and, 
both, thermal and rotation-induced radial growth [ 13 ]. 

Lubricant properties at the operating temperatures are 
evaluated once. Density is computed relative to a reference 
value at 16°C (60°F) 

P = p i6 * G (T c - 16) CD 

Kinematic viscosity (v) at atmospheric pressure is 
obtained from Walther's equation [ 14 ]. 

1Og 10 1Og 10 (V + * 6) = A* - B* log 10 (1.8T c + 492) (2) 

9 

SKF TECHNOLOGY SERVICES 

5KF : INDUS TPlfcS INC' 



4185C-I-R100 


AT81D006 


Dynamic viscosity is given by: 

n = (9.8 x 10 3 ) p\> (3) 

0 

The pressure viscosity coefficient is computed from a relation 
developed by Fresco [15]. 

(.^ 560 _) 

a = (3.3403 x 10' 8 ) [ C* ♦ D*log 10 v+ E* (log 1() v) 2 ] (4) 

Values for constant terms in equations (1) through (4), 
for two lubricants, are given in Table 1. 

2 . 2 Development of Mathematical Bearing Models 

The equilibrium field equation set, a global mathematical 
model, is formulated by summing the loads contributed by spe- 
cific local forces. Figure 7 shows those forces which have 
been considered in the analysis to act on a roller. Details 
of the mathematical models used to compute them follow. 

Roller to raceway contact elastic load is computed to 
recognize the influence of material properties, geometry, and 
component displacement. Given a set of values for the inde- 
pendent variables, component locations are established in 
space. Assuming Hertzian contact, and using the slicing 
technique [16]* relative positions of the rollers and rings are 
examined at each slice to determine if bodies contact or if 
free space exists between them. Relative position is estab- 
lished by using standard coordinate frame transformations 
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3 

[17]. ‘ Locations and magnitudes of interpenetration are 
recorded, with consideration for change in clearance, and 
well established equations [ 7 ] are used to determine maximum 
contact stress (a), contact half width (b) , elastic force (F) 
and moment (M) . 

Roller to raceway contact friction loads (fo and fj) are com- 
puted as functions of position and speed. Let £ be a vector locating 
a slice contact point, then the velocity of the roller surface 
at this point, relative to an observer moving at speed w 0 , is: 

v ■ a) r x p (S) 

The velocity of the raceway surface at the same point 
is : 

1 2 = (“* - £» 0 ) x E (6) 

where u>* = 0 for outer raceway contacts and w* = w for inner 
raceway contacts. Equations (5) and (6) enable computation 
of the sliding and entrainment velocities 


V = V - v 
— —2 —1 

(7) 

U = 1/2 (v + V ) 
— — 1 “2 

<- — -\ 

00 
v — / 


EHD film thickness is computed at each slice using the 
Loewenthal [18] model: 

3 

See Appendix A - "Coordinate Frame Definitions and 

Transformations" . 

11 
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h - K t (p') * (n |U|p 7E )' 62 (o/E') ^ 

o 0 — - 


(9) 


Lubricant inside tne concentrated contact will be exposed 
to a shear stress proportional to sliding speed. Assuming 
Newtonian fluid, linear velocity gradient, and exponentially 
varying fluid viscosity, the Allen [l9l traction model is used 


to obtain shear stress as 

t = n exp(acx) ( 2Z - lz ) 
0 


(10a) 


or 


T = f 

A 


when n exp(aa) ( V * z - Vlz ) v f A a 

0 — ^ A 

h 0 


and n exp(oto) ( 




1 z 


) > 


(10b) 


The magnitude of the contact friction force, |f | , is evaluated by 
integrating the shear stress, equation (10), over the contact 
area. Integration is performed using Simpson's Rule, assuming 
uniform velocity over the slice contact area and a semicylindrical 
pressure distribution. The friction force vector is directed 
along the sliding velocity vector 


= U! 


V 


lyi 


(in 


The moment generated by friction is given by 

m = p x f (12) 

Heat generation rate is taken as the thermal equivalent of 
mechanical work. 


q = f • V 
i — - 


(13) 


12 
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The EHD film starts in a gently narrowing inlet region 
where lubricant is entrained into the contact. Figure 7. 

Pressure builds up gradually in this hydrodynamic inlet, adding 
to rolling resistance. Computation of the net hydrodynamic 
force (h) and moment (£) is performed using the methods present- 
ed by Floberg [20] and Chiu [21]. Both cavitation and star- 
vation are considered. 

Roller end to flange load computation is similar to that 
made for the raceway. Independent variable values are used 
to establish, by coordinate transformation, the relative 
position of rollers with respect to the flange. Contact point 
location (p£) and penetration magnitude are computed in 

closed form [22] . 

The roller end to flange contact area, for sphere- ended 
rollers against an angled flange, can take on three possible 
shapes for a given load: 

(1) The contact area is Hertzian and elliptical 
in shape for relatively small end sphere 
radii over a limited, but practical, extent 
of roller to flange positions. 

(2) As end radius is increased, the contact area 
increases to a size bounded by the finite 
dimensions of the flange and roller end, pro- 
ducing high stresses at the boundaries. 

(3) When the sphere radius is large, the roller 
end may touch the flange at two non-Hertzian 
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contact points. 

Type (1) Hertzian point contact has been shown, by experiment 
[23], to be desirable because of its high load carrying 
capacity and ability to form a lubricant film. 

In the current work, Hertzian point contact at the flange 
is assumed. Contact maximum stress (o£) and dimensions 
(af»b£-) are evaluated as functions of geometry, material prop- 
erties and load [7]. 

After a program execution, an examination of predicted 
contact dimensions and contact point locations enable the 
determination if type (1) Hertzian contact exists. Subsequent 
executions can be made to select roller and/or flange geometries 
which provide this condition. 

The roller end to flange contact force (ff) is set so 
that the magnitude of its axial component, ff x , provides roller 
axial equilibrium. The f£ and f £ z components of f£ are com- 
puted as a function of sliding and entrainment velocities at 
the contact point |>£. Both EHD and inlet region HD force 
components are considered [19,20,21]. Computation of sliding 
and entrainment velocities is accomplished by using the point 
££ in place of p in equations (5) through (8). 

Moment (£f) is computed from the vector cross product 
Pf x f£, and heat generation rate is taken as the thermal 
equivalent of mechanical work 
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q 2 B if • { Cw - m 0 ) x p f - (w r x j> £ )} Ci4) 

Quasidynamic roller loading has been included to apply 
the equations of static equilibrium in the presence of accel- 
eration. Four specific terms are considered here. Roller 
centrifugal force (C) and the gyroscopic moment (G) vectors are evaluated 
as functions of roller position and speeds [24]. The two 
additional terms, included to reflect the inertial forces 
imposed during changes in roller orbital and rotational veloci- 
ties, are given by Kellstrom [25] as 


F 

-q 


_ Mn 
" FW 


i 

0) 

— o 


i + 1 

D AVG [ 1 2*0 ' 


i - 1 

I*. I ] 


k 


(15) 


M = MnD 2 |(o I 1 [ |a> l i + 1 - |(o | i ' 1 ] i (16) 

JTT ~ r ~ T 

where is the bearing's pitch diameter. 

Roller drag is included to account for forces generated 

when the lubricant is churned as it passes through the bearing 

cavity. The analytic description of roller drag presented 

here is obtained from equations describing external flow over 

cylinder [26] . 

2 

d = -9.8 x 10 3 (pc) D AVG C D 1 e D (w . w ) k (17) 

g -® 

The term pc in equation (17) is the effective density of 
the air-oil mixture in the bearing cavity, and must be chosen 
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to reflect bearing operating conditions. Correlation of experi- 
mental data with program predicted results [27,12], for a 40 mm 
jet lubricated spherical bearing, has shown c lies in the range 

.01 < c _< .02 

for shaft speeds to 5000 rpm. It is expected that this range o r 
values would be valid at higher speeds also. Note that bf. se 
spherical roller bearings are typically operated at low s\. , the 

drag force is also very low. The 40mm bore spherical roller bear- 
ing operating at 5000 RPM, for example, will experience a drag 
of . 3N (.07 lbs) per roller for C = 0.02. 

Roller to cage pocket contact load is a function of lubricant 
properties, speed, geometry and the position of the roller within 
the pocket. In the current work, the cage pocket geometry is 
taken as a rectangular cavity. 

It has been shown [3] that under typical operating conditions, 

skew moments (moment loading about the y axis, Figure 4) are 

IT 

imposed on the roller at both the outer and inner raceway contacts. 
Depending upon the magnitude of the skew moments, one of the 
following situations can occur: 

(1) The roller skews about the axis by some 
small amount. This is sufficient to balance 
the raceway induced moments with negligible 
assistance from the cage pocket or flange. 

(2) The roller continues to skew past 
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••he above condition until raceway in- 
duced moments are balanced by mutual con- 
tributions from the flange and cage 
pocket . 

(3) The bearing is manufactured without an 
integral flange. The roller continues 
to skew past the condition in (1), until 
equilibrium is achieved through balance 
of raceway and pocket induced moments. 

For he analytic assessment of pocket- induced skew moments, 
it is assumed that loads are sufficiently light to produce 
hydrodynamic contact, and that the pocket is fully flooded with 
lubricant. The slicing technique is then used to evaluate 
the forces and moments generated at the various stations along 
the skewed roller- to-pocket contact. 

In l28l, it was shown that roller position within the cage 
pocket could be written in terms of the roller’s speed and cage 
speed. Assuming the cage speed vector, oi , is coincident with 
the centerline of the .nner ring and its magnitude is the 
average of all roller orbital speeds, the j-th roller center 
to cage pocket center offset is 


(Z) 


j 


^AVG 77 
7n 


t'j 

Z 

1=2 



l-l 


+ 


2c 



1} + rz ) 1 


(18a) 


j = 2,3, 


,n 


u. 
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where 

(Z) 1 * ± 1/2 D avg sin ¥± sin 4>j + c 2 cos 4^ (18b) 

¥ * angular displacement of cage relative to its 

reference position 

= azimuthal angle of first roller 
c * cage radial displacement in y direction 

c z = cage radial displacement in z direction 

Note that equation (18) is applied to each row of rollers, 
where the upper sign refers to row 2 and the lower sign to row 1. 
The minimum hydrodynamic film thickness at the leading (+) and 
trailing (-) edge of each slice is given by (Figure 8) 

1 

(h+) 1 = - (Z) 1 + x sin (y s )* - r- cos (y s ) X (19) 

1 

(h-) 1 = + (Z ) 1 - x sin (y ) 1 - r- cos (Y.) 1 (20) 

0 4 5 a 5 

i 1 ^ 2 | • » • n ^ 

where r— is the radius of the slice loacated at x = x. 

Friction force acting on a particular slice, (R^)_, is com- 
puted as a function of the normal force. Equations ari used which 
describe the hydrodynamic lubrication of a rigid cylinder near a 
flat plate [29]. The normal force is obtained as a function of 
the film thickness, using a linear approximation of the hydrodyn- 
amic equations. 

Tne heat generation rate is taken as the thermal equivalent 
of mechanical work 
• 

(q,)_ - (R)_ (r_ IsU) C2i) 

X y X X 


SKF TECHNOLOGY SERVICES 

SKI INDUS IMS. INC ' 


18 



AT81D006 


The sum of the individual forces acting on each slice 
yields the total cage pocket force (R) and moment (S) . 

Cage rail/ring land torque and forces are computed from 
the hydrodynamic solution for short, self acting journal 
bearings [ 30] . 


T 



n r 
o 


(1-e 


3 

g 


w 


2 ) 1/2 



( 22 ) 


F = 

y 


n r w3 
o g g 


. W + CJ 


2,2 


(1-e ) 


(23) 


n r w 
o g g 


4 (1-e 2 ) 3/ 2 


iO + d) 


(24) 


Several bearing designs employ cages having multiple rails. 
Here, torque and radial loads are computed by specifying the 
rail land diameter and rail width of an equivalent single- 
railed cage. Dimensions are chosen so that the single rail 
equivalent cage and multiple rail existing cage produce 
identical torque and radial force at the same eccentricity 
values. The equivalent rail land radius is: 

r g = (ai/faVb)* 375 (25) 
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And the equivalent rail width is: 

w = (a 3 /b)- 125 (26) 

o 

where : 

l = number of rail' 
a = Z (rw^)» 
l = 1 

i = number of rails 
b = Z (r^w) 

£ = 1 4 

r = rail land radius of the 
£-th rail (M) 

w = the rail width of the 
l-th rail (M) 

Torque generated at the interface between split cage 
halves is computed by assuming a uniform pressure distribution 
over the interface and a constant coefficient of friction: 


T = 
— s 


+ 




(Fx)y 


(27) 


Where y is the coefficient of friction at the split, F x is the 
axial load supported at the split and A and B are radii defined 
in Figure 9. Torque computed by (27), is directed so that 
the slower cage piece is accelerated and the faster cage piece 
is decelerated. 

Heat generation rate at the split is given by: 
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^4 = I s • (“cf — C 2 ) C28) 

where and are speeds of each cage piece. 

Heat generation rate at the cage rail is given by 

q 5 - T . (^ - u) ( 29 ) 

2. 3 Formulation of Equation Set 

The field equation set is formulated by summing the loads 
acting on the rollers, cage and inner ring. 

Equilibrium equations for the i-th roller can be repre- 
sented in vector format by 2 equations 


n 


* [®k.m + fOil., * K.m * 1 * + 

m=l k=l ’ 9 9 


(fq ) 1 + (If)' + (D ) 1 = " 


( 30 ) 


and 


Z I 

m=l k= 1 





♦ (S)*] + (G) 1 + 


(Mq ) 1 = 0 ( 31 ) 

The cage maintains equilibrium through balance of the 

torques induced at the pockets, rails and split interface 

n 

D r n i i A 

(T r + T - AVC Z Z (- 1 J ) (R j) i = 0 ( 32 ) 

5 x 2 j=l i = 1 Z 
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Note that pocket loads are assumed to act at the pitch 
circle, and that if a split cage design is addressed, the 
additional term T g (equation (27)) is included; otherwise the 
term is zero. Summation of pocket and rail forces in the 
radial directions yields: 


i= 

?■ (k_)- sm <J>. + F cosg 
i =i L A 1 y 


sin8 = 0 


(33) 


i = n t 

£ (R )• cos<{>. + F sinB + F cosB = 0 (34) 


B is an angle measured CCW positive from the Y c axis to 
the point of closest approach between the cage rail and ring 
land. 


Three inner ring equations of force equilibrium are 
represented by ihe vector 


n t 

E 

i=l 


£ l®k,2 * <«k,2 * tWk.zl * ME) - o (35) 


where P - cage- land force vector for inner ring riding cage 
2 . 4 Computation of Planet Bearing Performance 

In the calculation of roller hearing performance it was 
assumed that the bearing rings are rigid 

and that elastic deformation occurs only at the Hertzian con- 
tact. In some special cases, such as planetary gear applica- 
tions, flexibility of the outer ring must he taken into account 
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in calculating the deflection at each roller location [6]. 

A typical, two row spherical planet bearing is shown in 
Figure 10. Note that the outer ring is integral with the 
gear to minimize component weight. 

Load is applied to the bearing at the (gear) pitch radius 
through diametrically opposed meshes. Figure 11. The bearing 
assembly maintains equilibrium by balancing the gear tooth 
loads with a reactive load at the carrier post. 

The roller to outer ring contact loads and gear mesh 
loads will tend to deform the normally circular shape of the 
outer ring. The deformed shape can be computed by superposing 
the influence of each individual component of applied load. 

The method is described in [6] , and requires that the gear 
mesh load, F, be represented by an equivalent set of loads 
applied at the outer ring neutral axis (Figure 12): 

R = F sin a 
T * F cos a 
M = F (Rp - R n ) cos a 

Assuming the deflections are elastic, the radial deforma- 
tion of the outer ring (measured as a displacement from the 
neutral axis) at the i-th roller location is 

A i * A i + A I " A i + A i (37) 
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where: 

M T R 

A A., A. - are the outer ring deformations due 
to moment, tangent and radial com- 
ponents of gear tooth load [29], 
measured positive inward at the 
i-th rolling element location. 

p 

A. - is the outer ring deformation due 

to rolling element loads, measured 
positive inward at the i-th rolling 
element location. 

The roller induced ring deformations can be written as 



C ij P j 


(38) 


4 

Where are influence coefficients for the outer ring and 
can be found in the literature [31] and Pj is the total roller 
to outer ring contact load of both rows at roller position j. 

Inserting (38) into (37) end rearranging terms, we obtain 
a set of i equations describing the deformed shape of the outer 
ring 


0 - ( Aj + 


.T * a* 
4 i * A i 


n 

E 


j = l 


C iJ P J - 


A i 


(39) 


Equation (39) is in a form which may be solved for 
deflections A^, using the Newton-Raphson scheme^. Note that 


The influence coefficients C- • are defined as the inward deformation of 
the ring at i due to an outwaTd unit load at j. 

^ When analyzing a spherical planet bearing, SPHERBEAN will formulate a set 
of deflection equations (39) and solve for along with the equations de- 
scribing inner ring, roller and cage equilibrium, (30) through (35). 
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the loads are non-linear functions of therefore equation 
(39) is non-linear. Also, in those kinematic cases where the 
carrier is moving, the additional centrifugal load on the 
rolling elements must be included in P^, Figure 13. 

The rolling element centrifugal load resulting from carrier 
rotation is computed by transforming the roller CG coordinates 
to a reference frame located along the sun gear center line. 

This enables computation of the centrifugal force vector, which 
is added to the centrifugal load due to roller orbital speed. 6 

2.5 Numerical Solution to the Formulated Bearing Analysis 

In contrast to the generality maintained in the bearing 
analysis formulation, specific solution procedures are required 
to address the specific problems arising during the iterative 
solution to equations (30) through (35) and (39) . 

An automated field equation set partitioning scheme was 
developed so that the convergence characteristics of several 
different subsets of equations and variables could be determined 
for a given load condition (i.e. pure thrust, pure radial or 
combined load) . A table containing 12 subsets which showed 
most rapid convergence was coded into the SPHERBEAN. Upon 
execution, SPHERBEAN examines the applied load to determine 

6 When analyzing planet bearings, "SPHERBEAN" will modify rolling element, 
cage and inner ring equilibrium equations to include proper centrifugal 
force. 
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direction and possible symmetries, and a corresponding pattern 
of equation subsets are selected from the table and sequentially 
brought to equilibrium using the Newton- Raphson method. 

Partial derivatives, required at each iteration of the 
Newton -Raphson method, are computed by fitting a second order 
polynomial through current and two adjacent function value 
points, and evaluating the derivative directly from the poly- 
nomial. Although more time consuming than Newton's forward 
difference operator, this procedure is employed to smooth the 
piecewise continuity of the field equation set and provide 
more reliable code performance. 

The numerical procedure is stopped when the RMS residual 
of all equations is less than a specified tolerance. Com- 
pletion of the iterative procedure permits the display of 

values for the vector of unknowns, bearing heat generation 

7 

rates and the computation of bearing fatigue life. The latter 
is computed according to standard Lundberg-Palmgren methods 
[32,33], and includes correction for the dependence of life on the 
film thickness to surface roughness ratio [34,35] . 

7 

Fatigue life computation details are given in Appendix B. 
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3.0 COMPUTATION OF SYSTEM STEADY STATE AND TIME TRANSIENT 

THERMAL PERFORMANCE 

Description to this point has addressed the bearing 
analysis at a given set of component temperatures. Thermal 
effects can be computed in either of two modes. 

The first mode details the steady state operating temper - 
atures of the bearing and its environment. SPHERBEAN will 
formulate the conservation of energy equations for an equivalent 
nodal model of the bearing and surrounding hardware. Energy 
equations contain heat sources which may be provided by the user 
to represent the heat generated at sources other than the 
spherical bearing, such as neighboring seals or gears. Energy 
equations will alsc contain the heat generated by the spherical 
roller bearing, (for example, see equations (13) and (14)). 

The Newton- Raphson solution method is used to solve the non- 

O 

linear conservation of energy equations for system temperatures. 
Resulting system temperatures are then used to compute bearing 
performance, including bearing heat generation rate. This 
iterative procedure, shown in Figure 14, is stopped when the 
difference between current and previous temperature predictions 
is less than a user specified tolerance (typically set at 2°C) . 


O 

The equation set is non-linear because it contains terms 
describing radiation and free convection (see Appendix C) . 
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The second mode can be used to detail time transient 
thermal performance of the bearing and its environment. Here, 
SPHERBEAN will formulate and solve a system of first order 
nonlinear differential equations. Typically, initial values 
are taken as the- solution to a steady state analysis. 

A detailed description of both steady state and time 
transient analysis can be found in Appendix C and in Volume II 
[11 1 of this work. 
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4.0 CONCLUSIONS AND RECOMMENDATIONS 

The material presented in the preceeding sections documents 
the creation of the state of the art spherical bearing analysis/ 
design tool SPHERBEAN. The formulations are general and the 
program architecture modular to permit easy maintenance and 


expansion. 

During the course of the work performed, a perspective 
was gained on the additional requirements for the accurate 
simulation of spherical roller bearing performance. 

Proper assessment of performance demands an investment 
in basic formulation and corresponding detailed experimental 
verification. It is specifically recommended that: 


1. The EHD traction and film thickness models be 
upgraded to include effects of surface micro- 
geometry. Because spherical roller bearings 
are typically operated at low speed, EHD film 
thickness at roller-to-raceway contacts will 
likely be several times less than the compo- 
site surface roughness in typical spherical 
bearing applications. 


2. The effects of lubricant application rate and 
distribution within the bearing cavity be in- 
vestigated. Thermal dimensional stability is 
sensitive to the assumptions concerning the 
lubricant distribution within the bearing 

cavity (C> and the sensitivity is accentuated 

9 

with increased speed. 

g 

See Example l in Volume II. 
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Since spherical roller bearings are being 
asked to follow the DN capabilities of ball, 
cylindirical and tapered roller bearings, 
realistic assessment of their high speed 
performance requires an accurate method of 
relating drag power loss to lubricant appli- 
cation rate. 

3. In light of the low film thickness to surface 
roughness ratios (A) that spherical roller 
bearings typically operate under, surface 
traction effects should be directly related 
to mater: al failure and thus bearing failure 
predictions . 

4. The kinematic description of inner ring dis- 
placement should be expanded to include a 
misalignment angle, and the performance of 
misaligned spherical roller bearings with out-^ 
er ring rotation be fully examined. 
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FIGURE 1 - SPHERICAL ROLLER BEARING 



FIGURE 2: SPHERICAL ROLLER BEARING GEOMETRY 
AND COORDINATE FRAMES 
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FIGURE 3: INNER RING DISPLACEMENT 
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FIGURE 4: ROLLING ELEMENT DISPLACEMENT 
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FIGURE 8: ROLLER TO CAGE POCKET INTERACTION 




FIGURE 


AT81D006 


CARRIER 

POST 


CARRIER 



4 - 


OUTPUT 


FIXED 

um / 


PLANET 

GEAR 

PITCH 

"RADIUS 



/777T7 

EXAMPLE OF KINEMATIC 
INVERSION WHERE CARRIER 
IS MOVING 


FIGURE 11: PLANETARY TRANSMISSION 
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FIGURE 15: (ITCH SPF.EO PLANF.T GEAR LOADING 




FIGURE 14: SIMPLIFIED FLOW CHART FOR SPHFRBFAN 
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APPENDIX A 

COORDINATE FRAME DEFINITIONS AND 
TRANSFORMATIONS 


u. 


SKF TECHNOLOGY SERVICES 


50 



-4J8SC .1 -R1Q0 


AT81D006 

APPENDIX A : 

COORDINATE FRAME DEFINITIONS AND TRANSFORMATIONS 

Consider a double row spherical roller bearing, Figure A-l > 
introduce a cartesian coordinate system - Y^ - Z n such that 
the X Q axis is coincident with the bearing's outer ring 
longitudinal center line. Position the Y q -Z 0 plane along this 
centerline such - that it contains the outer ring sphere origin, 
Figure A- 2. The Y q axis defines the angle 0=0° and the Z Q 
axis d> = 90°. Assume the X -Y -Z^ system to be fixed with 
respect to the outer ring, and the outer ring center of mass 
fixed in space. 

Introduce a second system, Xj-Yj-Zj, initially coincident 

with the X -Y -Z. system, but attached to the inner ring and 
o o 0 

free to move through space as the inner ring moves. Figure A- 3. 

A third coordinate system, X R -Y R -Z R , is introduced at each 

roller location so that the X R axis is along the roller center- 

line (rotation axis). The origin of the system is selected 

to lie on the line directed radially outward along the contact 

angle. Figures A-l and A-4. The Z R axis is tangent to the 

pitch circle. The X R -Y R -Z R coordinate system is fixed in space 

in a radial position so that the clearance between the roller 

and raceway, at the point of closest approach, is equal to 1/4 

the diametral clearance, Figure A-4. 

Introduce a fourth coordinate system, X_-Y -Z_, at each 

R E R 

roller location. This system is initially coincident with the 
X R -Y R -Z R system, but free to move through space as the roller 
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moves . 

RELATIONSHIPS BETWEEN COORDINATE SYSTEMS 

It is frequently necessary to express the components of a 
vector in at least two coordinate frames. This is particularly 
so in the definition of bearing geometry, where initial speci- 
fication is convenient in inertial coordinate frames. However, 
the bearing analysis frequently requires redefinition in 
frames which are in motion. 

The linear orthogonal transformation operator, [ 4>TS ] de- 
scribes the rotation portion of the transformat ’ Jn between 
"S" and "T" coordinate frames: 



*11 

*12 

*13 

[ 4>TS ] = 

*21 


*23 


*31 

*32 

*33 




— 

Six transformation operators will 

be < 

defined 


Cl) 


these all required transformations can be performed. They are 

1. Outer ring coordinate system (0) to the 
inner ring system (I). 

2. Outer ring coordinate system to the i-th 
roller inertial system (R) . 

3. i-th roller inertial system to the i-th roller 
moving system (TO . 

Transformations 4, 5, and 6 are the inverse of 1 , 2, and 3. 
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Consider the transformation from the outer to inner ring 
frames : 

I rn?v rnsv sinv -rn<:v sin\ 

T 


[$10J = 


cosy 2 

cosYy 

sm Yz 

-COSY z 

sinY 

-sin Yz 

COSYy 

COSY z 

sinY z 

sinY 

sinY y 


0 

COSYy 



( 2 ) 


where: Y z - rotation, or "misalignment an|le" of 

outer ring about the Z q axis . 

Yy - rotation of outer ring about the Y q 
axis . 

The inverse transformation, from the inner ring to outer 
ring frame, is given by the transpose of (2): 


[ *01 ] = 


cosy cosy 
z y 


smY, 


cosy smy 

z y 


•sinY z cosYy. sxnY, 


cosy 0 


smY z smYy cosy. 


(3) 


Similarlv, we obtain the following for the remaining four 
transformation operations: 


[ $R0 1 = 


cosa -sina cos<t>^ 

sina cosa cos<j>. 


0 -sin<J>^ 


-sina sin<|>- 
cosa sin^i 


cos4>. 


C4) 


1( ln 


the current version of SPHERBEAN, Y,, = Y_ = 0. 

y z 
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[ 4>0R ] = 


cosa 

-sina cos4>j 


-sina sintf^ 


sina 


cosa cos<J>^ 
cosa siniji- 


-sin<J>^ 
cos<J> • 


(5) 


[4RR ] = 


cosY t cosy s 

-sinY t cosy s 
siirr 


sinY, 


cosy, 


- cos y t siny s 
sinY t sinY s 
cosy_ 


(6) 


[$RR 


cosy, cosy 

L 5 

sinY„ 


-cosY t sinY s 


-sinY t cosy s sinY, 

< ~Y f 0 


->~nY t siny s cosy, 


(7) 


In equations (4) through (7) : 


• a, the contact angle, is defined as positive for 
rollers in j* 1 and negative for rollers in row 
1 , Figure A- 1 . 

* y is commonly referred to as the roller tilt 
'ingle, and represents the angular displacement of 

roller about the Z_ axis. 

R 

• y s is commonly referred to as the roller skew 
angle, and represents the angular displacement 
of the roller about the Y n axis. 
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A PPENDIX B 

FATIGUE LIFE COMPUTATION 
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INTRODUCTION 

Within SPHERBEAN, roller bearing fatigue life is calculated 
using Lundberg -Palmgren [32,33j methods. The value computed 
is then modified by multiplicative factors which account for 
material and lubrication effects. 


ROLLER BEARING RACEWAY LIFE 


The load distribution across a line contact is represented 
by a number of slices. The L^ 0 fatigue life of a given slice 

^cmk . 


is : 


L lOmk ' ( 


^emk 


) 


( 1 ) 


Here, Q cm ^ is defi ed as the load for which a slice will have 

90 percent assurance of surviving 1 million revolutions. 

Letting m refer to a raceway, k to a slice, where n is the 

index of the last slice, the explicit form of Q ^ is given by [8] 

49500 X { D ( 1 - Y m ) }1 * 074 * m °' 778 (Y m )°’ 222 f 2) 

^cmk fir.-i "+ -,1 0. 25 


tN(l - Y m )} 


where: 

N: is the number of rollers per row 
X: is the capacity reduct >n factor 

X^ = .8 when k = 2,3,.. , (n-1) 
X^ = .53 when k = 1 or k = n 
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Note that the upper sign is used for the outer race, the lower 
sign for the inner race. 

Q em k i s the equivalent load for the slice: 




Q m kj is the individual contact load on the k-th slice and 
e = 4.0 or e - 4.5, depending respectively upon whether the 
applied load rotates or is stationary with respect to the race- 
way in question. 

L 10 life of a raceway is given by 

n -e -1/e 

L 10m = a 2 a 3 a 4[ k f 1 ( L lO mk ) } 

where e is the Weibull slope exponent, here taken to be 9/8, a: d 

a2 is a life improvement factor to account for 
improved materials (User input to program, see 
[ 35 ] factors ”D”and "E.") 

a^ is a life factor to account for film thickness 
to surface roughness ratio (computed within 
SPHERBEAN [ 34l ; factor "F "in [35] ). 

is a factor which accounts for materials having 
a modulus of elasticity other than that of basic 
steel (computed within SPHERBEAN). 

BEARING LIFE 

Ljq life of a single row bearing, considering both raceways, 
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is : 


2 

E 

m= 1 


■10 ■ < J, < L 10m’ 1 


e, -1/e 


(5) 


If the bearing contains two rows of rollers, each row can be 
treated as a separate bearing and the lives summed to yield 
the double row bearing life 

-1/e 


L 10 " { 


10 


ROW 1 


) 


(L 


10 


ROW 2 


)' 6 } 


( 6 ) 
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APPENDIX C: THERMAL MODEL 
TEMPERATURE CALCULATIONS 

After each calculation of bearing generated heat rates, 
either steady state or time transient temperature analysis may 
be performed. The computations are terminated in the following 
manner: 

1. The steady state case terminates when each 
system node temperature is within A Centi- 
grade of its previously predicted value. 

The value for A is specified by the user 
(typically 2°C). 

2. The transient calculation terminates when 
the user specified transient time interval 
is reached or when one of the system temp- 
erature nodes exceeds 600°C (1112°F). 

STEADY STATE TEMPERATURE MAP 

The physical structure is considered to be divided into 
a number of elements represented by nodes. Heat fl to node 
i from surrounding nodes j, plus the heat generated at. node i, 
must equal zero to satisfy the definition of steady state 
conditions . 

After each calculation of bearing generated heat, which 
results from a solution of the bearing portion of the program, 
a set of system temperatures is determined which satisfy: 

q i = q Qi + q gi = o for all nodes i (1) 
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where q Q ^ is the heat flow from all neighboring nodes to node i 
q • is the heat generated at node i. Values are 

o 

calculated to represent heat created by bear- 
ing friction. 

The resulting set of field equations is solved with a 
modified Newton- Raphson method which successfully terminates 


when 


n CEQ) ? 1/2 < . 

£ [__L. ] 5 


( 2 ) 


where, 

n = number of field equations 
EQ^ = residue of the i-th field equation 
6 = .01 

TRANSIENT TEMPERATURES 


The net heat q^ transferred to the i-th node heats the 


element , i . e. : 



(3) 


where 

P = density 
C = specific heat 
V = volume of the element 
t = temperature 
T = time 

The temperatures, t at the time of initiation T=T g are 
assumed > be known, that is 
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t oi = *i ( V 1 = (*) 

The problem of calculating the transient temperature dis- 
tribution in a bearing configuration thus becomes a problem 
of solving a system of nonlinear differential equations of 
the first order with prescribed initial conditions. The equa- 
tions are nonlinear since they contain radiation terms and 
free convection, which are nonlinear with temperature as will 
be shown later. The simplest and most economical way to arrive 
at a solution is to calculate the rate of temperature increase 
at the time T = T^ from equation (3) and then compute the 
temperatures at time T^ + AT from 


dt, 


t k+l t k + 


dT 


AT " t k + p rr 


AT 


(5) 


CALCULATION OF TRANSIENT TIME STEP . 

If the time step AT used is chosen too large, the tempera- 
tures will oscillate; if it is chosen too small, the calculation 
will be costly. It is therefore desirable to choose the 
largest possible time step that does not give an oscillating 
solution [36, 37] *. 


dt- 


dt • 


i ,k+l >n -! — i o 

— 9 0 i = 1, 2, ..., n 


i,k 


( 6 ) 


u. 
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If this derivative were negative, the implication would 
be that the local temperature at node i has a negative effect 
on its future value. This would imply that the hotter a region 
is now, the colder it will be after an equal time interval. An 
oscillating solution would result. 

Differentiating (5) for node i and combining with (6), the 
time step size condition is 

dt- AT i dq i 

" 1 + "p.C V. 3tT >jr 0 (7) 

i Pi i 


The derivative dq^/dt^ is approximated using the forward 
difference operator 


d qi 

(Tt~ 


Ct ± + At A ) - qiftj) 


( 8 ) 


At i 


The values AT- which satisfy the equality in equation (7) 
are calculated. The array is searched, and a value of AT , 
rounded off to one significant digit smaller than the smallest 
of the ATi obtained is used. 

CALCULATION OF HEAT TRANSFER RATE 

Heat transfer mechanisms which occur in a bearing appli- 
cation are: 

• Conduction between inner ring and shaft 
and between outer ring and housing 

« Convection between the surface of the 
housing and the surrounding air. 
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• Radiation between the shaft and the housing. 

, Forced convection between the bearing and cir- 
culating oil. 

All the above mentioned modes of heat transfer are consid- 
ered in calculations of the heat balance at each given node. 

GENERATED HEAT 

A heat source may exist at node i. The quantity representing 
the source magnitude must be added to the net heat flowing from 
neighboring nodes. 

When the heat source is other than a spherical roller 
bearing, it may be considered to produce known amounts of 
power, in which case constant numbers are entered as input to 
the program (see Example 1 in [11]). 

CONDUCTION 

The heat flow q • . which is transferred by conduction 

Cl 9 J 

from node i to node i, is: 


'cl 



( VV 


where A = the thermal conductivity of the medium 
l - length between i and j 

FREE CONVECTION 

Free convection between a solid medium and air, the heat 
flow q . • transferred between nodes i and j can be calculated 

V 1 y J 

frrn the equation 
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VI, J 


= a A 
v 


iv^r 


SIGN (tj-tj) 


( 10 ) 


where 

a y = the film coefficient of heat transfer by 
free convection 

A = the surface area of thermal contact between 
the media 

d = is an eirwnent, usually = 1.25, but any 
value can be specified as imput to the 
program 


SIGN = 


1 if ‘i > ‘j 


1-1 if t < t, 
i J 


The value of can be calculated for various cases [ 36, 38] . 
FORCED CONVECTION 

Heat flow q - - transferred by forced convection can be 

Wl y J 

obtained from the following equation. 


q wi,j “w A ^i ’ 


(ID 


where a is the film coefficient of heat transfer during 
forced convection. This value is dependent on 
the actual shape, the surface condition of the 
body, the difference in speed, as well u j the 
properties of the liquid or gag [38]. 

In most cases, it is possible to calculate the coefficient 


of forced convection from a general relationship of the form, 

N u - a R V (12) 
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where a, b, and c are constants obtained from handbooks, [39 ] , 

R and P are dimensionless numbers defined by 
e r 

N u = Nusselt number = a w L / X 
L = characteristic length 
X * conductivity of the fluid 
R g = Reynold's number = ULP/n 
U = characteristic speed 
p = density of the fluid 
n * dynamic viscosity of the fluid 
P r = Prandtl's number = nC p /X 
Cp = specific heat 

SPHERBEAN can accept a specified constant value for the 
coefficient of convection, or, at the user's option, the 
coefficient can be calculated internally by the program. If 
the calculation option is exercised, input can be given in 
one of three ways: 

Constant Viscosity 

1. Values of the parameters in Equation (12) are 
given as input and a constant value of a w is 
calculated by the program. 

Temperature Dependent Viscosity 

2. The coefficient a for turbulent flow and heating 

w 

of petroleum oils is given by 
a w = kg . { n(t)> k 10 

where kg and k^g are given as input together 
with viscosity at two different temperatures. 
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3. Values of the parameters in Equation (12) 
are given as input. Viscosity is given at 
two different temperatures. 

RADIATION 

If two flat, parallel surfaces of same surface area A, are 
placed close together, the heat transferred by radiation between 
nodes i and j representing those bodies, will be, 

q Ri,j = eCTA [(t i + 273)4 " Ct j + 273 ) 4 ) (I 4 ) 

here, e is the surface emissivity, and o is the Stefan-Boltz- 
mann radiation constant. 

Heat transfer by radiation under other conditions can also 
be calculated [36 ,38 ] .The following equation, for instance, 
appl ies between two concentric cylindrical surfaces: 

eaA i [ (t i + 273) 4 - (t- + 273) 4 ] 

q — 

’■> 1 + (1-e) (A i /A e ) 

where A^ is the area of the inner cylindrical surface 
A g is the area of the outer cylindrical surface 

FLUID FLOW 

Between nodes established in fluids, heat is transferred 
by transport of the fluid itself and the heat it contains. 

Figure C-l, on the following page, shows nodes i and j 
at the midpoints of consecutive segments established in a stream 
of flowing fluid. 
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FIGURE C-l: FLUID HEAT NODES 


The heat flow q . . through the boundary between nodes 

U1 y J 

i and j can be calculated as the sum of the heat flow q fi 
through the middle of the element i, and half the heat flow 
q Q ^ transferred to node i by other means, eg., convection. 
The heat carried by mass flow is, 


q fi * p i C Pl V i *i = "i'i 


06) 


where = the volume flow rate through node i. 

The heat input to node i is the sum of the heat generated 
at node i (if any) and the sum over all other nodes of the 
heat transferred to node 1 by conduction, radiation, free and 
forced convection. 


q oi = q G,i + J ^ q ci,j + q vi , j + q wi,j + q Ri,j 
The heat flow between the nodes of Figure C-l is, then, 


q D < O (17) 
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q ui.j * «fi * v /2 (18) 

If the flow is dividing between node i and j, as in 
Figure C-2, then the heat flow is calculated from 


V,j ’ K ij (q 


fi ♦ <W 2 ) 


(19) 


where = the proportion of the flow at i going to node j , 

04 K- . * 1. 



TOTAL HEAT TRANSFERRED 


The net heat flow rate to node i can be expressed as. 


( 20 ) 


The summation should include all nodes j, which interact 
with i. Unknown temperatures as well as those specified as 
known should be included. 
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conduction through a bearing 

The conduction between two nodes is governed by the thermal 
conductivity parameter X. The value of A is specified at 
input. 

An exception occurs when one of the nodes represents a 
bearing ring and the other a set of rolling elements. Here, 
the conduction is calculated separately by the program using 
the principles described below. 

THERMAL RESISTANCE 

It is assumed that the rolling speeds of the rolling 
elements are so high that the bulk temperature cf the rolling 
elements is the same at both the inner and outer races, except 
in a volume close to the surface. The resistance to heat low 
can then be calculated as the sum of the resistance across 
the surface and the resistance of the material close to the 
surface. 

The resistance is defined implicitly by 

At * R • q (21) 

where At is the temperature difference 
and q is the heat flow 

The resistance due to conduction through the EHD film 
is calculated as 

Rj * (h/A) . A (22) 
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where h is taken to be the calculated plateau film 
thickness 

A is the Hertzian contact area at the specific 
rolling element-ring contact under consideration. 

X is the conductivity of the oil. 

The geometry is shown in Figure C-3 (a). Asperity conduc- 
tion is not considered. 

So far, a constant temperature difference between the sur- 
faces has been assumed. Pul during the time period of contact, 
the difference will dec e because of the finite thermal 


diffusivity of the materia", near the surface. Figure C-3(b). 

To points at a distance from tne surface, this phenomenon 
will have the same * ' feet as an additional resistance ^ act ~ 
ing in series with 

This resistance was estimated [40 ] as. 


re. 1 


, ** %l/2 

( iF7r _) 


where i = contact length, or in the case of ;.n 
re 

elliptical contact area, 6.8 times 
the major axis. 


X = heat conductivity 


ip = thermal diffusivity = X/ ( p. C^) 
p = density 


Cp = specific heat 
b = half contact width 


V = rolling speed 

7; 
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The resultant resistance is 


"res ’ n l * «2 


(24) 


There is one such resistance at each rolling element. 

Thev all act in parallel. The resultant resistance, ^ reS » is 
thus obtained from 


1 - " 1 
fi res i=l ^res,i 


(25) 
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[a,i Schenatic concentrated Contact 




direction of rolling 


(b) Temperature Distribution at Rolling, 
Concentrate' 4 Contact Surfaces 


FIGURE C-3 : CONTACT GEOMETRY AND TEMPERATURES 
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